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The elastic electron-proton scattering data are analysed in order to determine proton charge 
and magnetic rms radii, j*e and ru- Along with the usual statistical error, we try to estimate a 
systematic error in the radii, caused by the inadequacy of particular form factor parameterization 
employed. The range of data to use in the analysis is chosen so as to minimize the total (statistical 
+ systematic) error. We obtain r E = 0.912 ± 0.009 (stat) ± 0.007 (syst) fm and r M = 0.876 ± 
0.010 (stat) ± 0.016 (syst) fm. The cross-section data were corrected for two-photon exchange. We 
found that without such corrections obtained rg and tm are somewhat smaller while the quality of 
fit is worse. 

I. INTRODUCTION 

Charge and magnetic rms radii are important low-energy characteristics of proton, showing its "size" . They are 
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where Ge and Gm are electric and magnetic form factors (FFs) and fi = Gm(0) is proton magnetic moment. In the 
nonrelativistic approximation the radii can be expressed via moments of charge and magnetization densities: 
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Dh" r| = Jp E (r)r 2 dV, r 2 M = / p M (r)r 2 dV (2) 
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The knowledge of proton radii is necessary for the calculation of other low-energy observables with taking finite proton 
size into account. The examples are Lamb shift and hyperfine splitting in hydrogen, where the experimental errors 
are several orders less than the proton finite size corrections [l| . 

Proton radii are usually determined from the elastic electron-proton scattering data. One chooses some FF param- 
t-H ' eterization and then performs a fit to experimentally measured FFs or cross-sections. The radii are then obtained by 
ON i Eq.([T]). Nevertheless, this procedure is somewhat ambiguous. At first, there are many ways to parameterize FFs. The 
most fundamental way is to construct a parameterization which has proper analytic structure (poles and cuts in the 
time-like region, asymptotic behaviour, etc.) @, H|. The drawbacks of this method are very complicated expressions 
t-H ■ for FFs, need for a model input, and need to analyze proton and neutron FFs simultaneously and for all Q 2 . If we 
are interesting just in proton radii, it is simpler to restrict oneself to electron-proton scattering at small Q 2 , as we 
do in the present work. The FF parameterization is not constrained anyhow; it may be chosen for the reasons of 
convenience and accuracy of results. The examples are linear, quadratic [J], continued- fraction [5[ parameterization, 
etc. 

Our work is most similar to Ref. @. As it was done in Ref. Q, we parameterize electric and magnetic FFs simul- 
taneously and fit the cross-section data directly, instead of performing Rosenbluth separations and fitting obtained 
FFs. However, we use another, more intuitively understandable, FF parameterization, normalize the data from each 
experiment independently (the normalization factors become fit parameters) , and evaluate both electric and magnetic 
radii. 

Besides the cross-section data, today there is a number of data on polarization observables. However we do not 
use the latter in our analysis, because the polarization observables do not depend on Ge and Gm separately, but 
depend on the Ge/Gm ratio only. The data on Ge/Gm can constrain a small difference r 2 E — r M , but not or Tm 
itself. Additionally, there are about 10 polarization data points for Q 2 < 1 GeV 2 , to be compared with hundreds of 
cross-section points available. Therefore the effect of polarization data inclusion would be extremely weak. 

The distinctive feature of our work is the new approach to the error estimation. Usually, two sources of the error 
are considered: 1) the statistical errors in the individual data points, and 2) the overall normalization uncertainty 
for each experiment. But there is third error source: the inadequacy of chosen FF parameterization. Indeed, if we 
try to use, for instance, linear parameterization for < Q 2 < 1 GeV 2 , we will never obtain accurate results, since 
the FFs are very different from a linear function on this interval. On the other hand, choosing too small Q 2 range 
increases first and second type errors, since less data points will be used and FF value will be too close to 1. The 
latter circumstance was also pointed out in Ref. [5], but to our opinion, author did not treat the third type of error 
properly. In particular, author studies the scatter of te, associated with the choice of fit type (i.e. number of terms in 
the parameterization and range of data to use), and counts it as a part of the statistical error. Instead, we assert that 
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fit type should be chosen so as to minimize the total error, which consists of three above-mentioned parts. Therefore 
the errors in Ref. [j| turn out to be overestimated. The detailed discussion of the error estimation is given in Sec. IIIII 
The importance of two-photon exchange (TPE) corrections for the proton radius extraction was several times 
stressed in the literature. We apply TPE corrections to the cross-sections before fitting and check their influence on 
the results. 



II. FITTING PROCEDURE 



The following x 2 function is minimized: 
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where e enumerates experiments, i — data points within each experiment, <7g* p and da e ^ are experimental cross- 
sections and their point-to-point errors, N e are normalization factors (to be determined by fitting), dN e are quoted 
normalization uncertainties and cr th is theoretical cross-section: 

a th = eG 2 E (Q 2 ) + Q 2 G 2 M (Q 2 )/4M 2 + <J^(Q\e). (4) 

Here a 2l (Q 2 ,e) is TPE correction, calculated according to Ref. [6J. Since we are interesting in proton radii, which 
are related to FF behaviour near Q 2 — 0, it is enough to take into account only the elastic intermediate state in 
the calculations of TPE. The inclusion of TPE corrections is important, since they not only do not vanish but grow 
rapidly as Q 2 — > 0, what may noticeably affect the radius extraction [7|. 
The FFs were parameterized as 

n n — 1 

G B (g 2 ) = (l-^o) 2 E a ^' GM(Q 2 ) = (l-£/£o) 2 ]>>£ fc ' (5) 

k=0 k=0 

where £ = Q 2 /(l + Q 2 /£ ) with £ = 0.71 GeV 2 . 

The multiplier in front of the sum is equal to (1 + Q 2 /£o) -2 and represents dipole falling of FFs, thus the remaining 
polynomial should be approximately constant over the full Q 2 range; this reduces statistical errors in its coefficients. 
The polynomial in £ is better than the polynomial in Q 2 , since £ remains finite at Q 2 — > oo, thus the behaviour 
of polynomial at higher Q 2 is well-controlled, whereas the polynomial in Q 2 inevitably grows with Q 2 . We use the 
polynomial of degree n for parameterization of Ge and n — 1 for Gm , because the contribution of Gm to the cross- 
section is damped by Q 2 /AM 2 [see Eq. (j4j], and thus the Gm form factor is less significant at low Q 2 . The choice of 
the parameter n is discussed in Sec. HVl 

The proper FF normalization is ensured by keeping a = Gb(0) = 1 and b = Gj/(0) = 2.793 fixed. So the free 
fit parameters are coefficients a\, a„; bi, 6 n — l an d normalization factors N e . After the fitting, te and Tm are 
calculated via Eq.([T]). 



III. ERROR ESTIMATION 



The statistical error can be computed in the straightforward way (see exact formulae below), and includes both 
the error caused by point-tojDoint uncertainties and the error caused by an uncertainty in the normalization of each 
experiment. Some authors Q identify the last one as a systematic error, but we do not agree with them. While the 
normalization errors are systematic with respect to a single experiment data, they become statistical with respect 
to the proton radius, which is extracted from a number of experiments, every one having its own normalization. 
Moreover, if an experiment is re-done from scratch, the normalization coefficient will be different each time, similarly 
to statistical errors. 

Instead, we recognize another source of really systematic error in the radius. Suppose for simplicity that FF is 
measured directly and we fit it with, say, quadratic polynomial. If we knew a priori (by some theoretical arguments) 
that FF in any case must have this definite form 1 , the only thing to do was to extract the coefficients of the polynomial, 



For example, if instead of FFs we were studying time dependence of radioactive decay rate, we knew that it must have exponential 
form. So in this case we should fit the data with exponential, not polynomial or the like. 
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and the error in the FF would be the consequence of the errors in these coefficients. But actually we do not know 
anything about FFs a priori (except the normalization at Q 2 — 0). Therefore even if (accidentally) the experimentally 
measured FF values coincide with the true ones, the FF obtained by the fitting procedure will deviate from the true 
one, simply because the assumed parameterization cannot reproduce true FF well enough. The resulting difference 
has all properties of a systematic error: it will not change if the experiments are re-done and it does not depend on 
the precision of the data. 

From the mathematical point of view, the radius rm, yielded by the fitting procedure, is a function of experimentally 
measured cross-sections cr^ p , 

rst = r fit (cr°7) (6) 

The measured cross-sections deviate from the true ones, thus 

*lJ = (l + AN e y™ + A* eti (7) 

where AiV e and Aa ei are unknown random quantities. The squares of quoted statistical and normalization errors, 
da Et i and dN e , used in Eq.Q, are variances of these quantities: 



dai ti = Aal v dN* = AN* (8) 

(the bar denotes mathematical expectation). Since Act and AN are small enough, we may expand rat into Taylor 
series about tr* 1 '" : 

rat =r + Y, CeA^eTlT + Aa e ,i] + ]T d e , i;e , y[AN e al™ + Aa^AN^a'jy + Aa c7 ] + 0(Aa, AN) 3 , (9) 
where 

~ ~ /_truc\ <9ffit / truc\ j \ 3 rfit / true\ {~\(\\ 

ro=rfit((J ), c e , a = d^cxp {o- ), a e ,»;e',i' = 2 da c *? da *? 
The quantity rp in general is not equal to true radius r truc ; the difference is what we call systematic error: 



dr syst =r - r truc = r fit - r true (11) 
The total error of the radius is naturally defined by 



dr 2 = (r fit - r true ) 2 (12) 
Using Eq.Q and assuming, as usual, no correlation between different data points and experiments, 

Ao- eii Aov,i' = S^S^dali, AN e AN e , =5 ee ,dN 2 , (13) 

we easily obtain 

dr 2 = (r - r truc ) 2 + ^ [(r - r tl . uc )d e ^, + c^c^} [o^offidN? + S^da 2 ^} + 0(dN, da) 4 . (14) 

If c?r syst is reasonably small, then the first term in the square bracket can be neglected; also cr truo can be safely replaced 
by cr cxp , thus we obtain 

dr 2 = dr 2 yst +dr 2 tat , (15) 

where 

dr 2 tat = ]T ., Kfa^dN 2 + 5 u ,daU (16) 



is the statistical error. 

The systematic error is not so easy to estimate as the statistic one. It follows from Eq.©, that the systematic error 
is the error of the radius that would be obtained if tr cxp coincide with cr truo . Thus to estimate it we need to know true 
FFs, but this is exactly what is unknown. The best t hing we can do is to use pseudo-data. First, we choose some 
realistic FF parameterization (from one of Refs. [2l.l5l.l8l-ll0| or dipole fit). After that rt ruc is known. Then we generate 
pseudo-data at e and Q 2 of experimental data, assign them the same uncertainties da e< i, but no errors {Aa e ^ = 0, 
A7V e = 0). Then the fitting procedure, described in Sec.|TIl is applied and the radius rgt is obtained. The systematic 
error is the difference rg t — r true , since the statistical errors are absent by construction (of course its value depends on 
the parameterization chosen). The maximal absolute value among different parameterizations is taken as an estimate 
of the systematic error. 
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IV. CHOICE OF FIT TYPE 



Before we perform the fitting of actual data, we must choose the polynomial degree n and the higher bound for 
momentum transfer Qj 2 nax (a data point will be included in the analysis if it has Q 2 < Q^ax)- We will show that both 
n and Q^ ax can be chosen to minimize the total uncertainty of the radius. 

Suppose that the polynomial degree n is fixed and Q^ax i s varied. The statistical error will decrease as Q, 2 nax 
increases, simply because more data points are used. On contrary, the systematic error will grow: the polynomial of 
some fixed degree may be a reasonable approximation of FF on a small interval of Q 2 , but becomes less adequate on 
a larger interval, because the actual FF is not a polynomial at all. 

Thus we suggest that for each n there is some optimal Q^fn) with which the total error is minimal. Indeed, this 
was confirmed by fitting of pseudo-data. As an example, the dependence of the error on Q 2 lax is shown in Fig. [I] for 
n = 2. The curves are not smooth because of discrete nature of the data. The optimal Q 2 nax (H) values are summarized 
in Table [I] These optimal values are different for electric and magnetic radii. This means that one does not obtain 
the best estimates for both te and tm from a single fit; it is necessary to perform two fits with different Q max — one 
for te (upper half of the table) and another for tm (lower half). 
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FIG. 1: Estimated errors in electric radius vs. Q max for n = 2, statistic (dash-dotted), systematic (dashed) and total (solid). 
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O 2 

max 


points 




dr Byst 


dr 


r 


X 2 /d.o.f. 


r 


X 2 /d.o.f. 




2 


0.10 


159 


0.0096 


0.0068 


0.0117 


0.895 


0.52 


0.885 


0.56 




3 


0.30 


261 


0.0088 


0.0066 


0.0110 


0.912 


0.70 


0.904 


0.74 




4 


0.98 


358 


0.0097 


0.0054 


0.0111 


0.914 


0.75 


0.905 


0.79 




2 


0.07 


126 


0.0175 


0.0248 


0.0304 


0.836 


0.50 


0.803 


0.52 


tm 


3 


0.19 


226 


0.0129 


0.0158 


0.0204 


0.885 


0.68 


0.859 


0.73 




4 


0.98 


358 


0.0100 


0.0158 


0.0186 


0.876 


0.75 


0.857 


0.79 



TABLE I: Fit parameters and results. Qmax values are in GeV 2 , r and dr — in fm. 



V. RESULTS 



Finally, the fitting procedure was applied to the actual experimental data, with n — 2. A and with optimal Q^^x 
for each n. The data were taken from Refs. [TT14T9| . The resulting radii and y 2 per degree of freedom are shown in 
the "Results" section of Table Q] Choosing different n we obtain somewhat different radii, but they agree with each 
other within their estimated uncertainties. The most precise values are obtained: for the electric radius with n = 3, 
te = 0.912 ± 0.011 fm, for the magnetic radius with n = 4, tm — 0.876 ± 0.019 fm. The magnetic radius is smaller 
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TE 


TM 


Ref. 


Method used 


0.847 


0.836 


P. Mergell et al. \2] 


eN - 


> eN, dispersion analysis 


0.822..0.835 


0.843..0.852 


M.A. Belushkin et al. [3] 


eN - 


> eN, dispersion analysis, pQCD approach 


0.840..0.852 


0.849..0.859 


M.A. Belushkin et al. [3] 


eN - 


■> eN, dispersion analysis, SC approach 


— 


0.778(29) 


A.V. Volotka et al. [20] 


HFS 


in hydrogen 


0.83(5) 




I. Eschrich et al. [21] 


pe — > 


ep 


0.880(15) 




R. Rosenfelder [22] 


ep — > 


ep, Coulomb corr. 


0.895(18) 




I. Sick [5] 


ep — > 


ep, Coulomb corr., continued fraction param. 


0.912(11) 


0.876(19) 


This paper 







TABLE II: Recent determinations of proton radii 



than the electric one, but has larger uncertainty (because the magnetic FF weakly affects low-Q 2 cross-section). In 
Table HI] our results are compared to other recent determinations of proton radii. 

Our best value for electric radius is larger than all others. However, it is still in agreement with the one obtained in 
Ref. [5] using continued fraction parameterization for FFs. Though our value is somewhat larger, their difference lies 
within estimated errors. Thus the difference does not imply that either value is incorrect, but gives an impression of 
the precision to which we know Te- Nevertheless, there is one interesting implication. The theoretical value of classic 
2s — 2p Lamb shift with proton size correction can be written as [l[: 

Ai/ = 1057687(2) + 196r|, (17) 

where Az/ is given in kHz and rg — in fm. With our te value this equation yields A^ = 1 057 850(2)(4) kHz, in 
agreement with the most precise Lamb shift measurement result, 1057 857.6(2.1) kHz [ll l23l|. 

To study the influence of TPE on the obtained results, we repeat the fit without TPE corrections to data [putting 
<j 27 = in Eq.Q]. The results are shown in the last two columns of Table |U We see that x 2 1S slightly bigger in 
this case, that is, TPE-uncorrected data are fitted worse. Both electric and magnetic radius increase after taking 
TPE into account, te by ~ 0.01 fm and Tm by ~ 0.03 fm. The change is comparable with the error for te, but is 
more significant for rjyf. This result refutes our claim made in Ref. Q, where we suggest, basing on Q 2 -dependence of 
TPE amplitudes, that the inclusion of TPE should lessen obtained radii. This is because the data points from typical 
experiment do not lie at fixed e, but lie at fixed energy. The interplay of s- and Q 2 -dependence leads to the increase 
of the radii. 

Though our work was aimed at the extraction of te and r&f, as a by-product we also obtain FF parameterizations 
in the low-Q 2 region. It may be interesting to compare them to other existing parameterizations, in particular, to the 
latest "global fit" [lCf. In Fig. [5] we plot the global Ge parameterization from Ref. [10( and three parameterizations 
arising from three "best fits" in Table |U Each of them is valid for Q 2 less than corresponding Q 2 ^^. All three 
curves behave similarly and differently from the fit of Ref. [Io[ . After the examination of data sets used in Ref. [l(| , 
we suppose that the difference is due to the omission of data from Ref. [l4j], included in the present analysis. The 




FIG. 2: Electric FF parameterizations, normalized by dipole FIG. 3: The cross-section at e = 1. Odipoic is the cross- 

FF Gdipoic = (1 + Q 2 /0.71)~ 2 . section calculated assuming dipole FFs. 
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- Fit, n=3 




Q 2 , GeV 2 
FIG. 4: hGe/Gm ratio. 



kinematics of this experiment is: e ~ 1 and Q 2 = 0.05 to 0.3 GeV 2 , thus the data impose important constraints on 
low-Q 2 Ge behaviour. The parameterization from Ref. [10| totally fails to describe these data, while in our analysis 
they were fitted reasonably well (see Fig. [3]). In Fig. [4] we plot ^jlGe/Gm ratio, derived from our FF parameterizations. 
This ratio is unity at Q 2 = and is known to fall significantly at Q 2 > 1 GeV 2 . However its precise behaviour at low 
Q 2 it not well established. Our results suggest that this ratio crosses unity at Q 2 ss 0.3 GeV 2 and slightly grows at 
lower Q 2 , returning back to 1 at Q 2 < 0.1 GeV 2 . It is interesting to note that similar behaviour may be favoured by 
recent polarization transfer measurements [24| . 
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